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Nuclear superfludity in exotic nuclei close to the drip lines and in the inner 
crust matter of neutron stars have common features which can be treated with 
the same theoretical tools. In the first part of my lecture I discuss how two such 
tools, namely the HFB approach and the linear response theory can be used to 
describe the pairing correlations in weakly bound nuclei, in which the unbound 
part of the energy spectrum becomes important. Then, using the same models, 
I shall discuss how the nuclear superfluidity can affect the thermal properties 
of the inner crust of neutron stars. 



1. Introduction 

The basic features of nuclear superfluidity are the same in finite nuclei 
and in infinite Fermi systems such as neutron stars. Yet, in atomic nuclei 
the pairing correlations have special features related to the finite size of 
the system. The way how the finite size affects the pairing correlations 
depends on the position of the chemical potential. If the chemical potential 
is deeply bound, like in stable and heavy nuclei, the finite size influences the 
pairing correlations mainly through the shell structure induced by the spin- 
orbit interaction. The situation becomes more complex in weakly bound 
nuclei close to the drip lines, where the chemical potential approaches the 
continuum threshold. In this case the inhomogeneity of the pairing field 
can produce a strong coupling between the bound and the unbound part of 
the single-particle spectrum. This is an issue which will be discussed in the 
first part of my lecture. More precisely, I shall discuss how the continuum 
coupling and the pairing correlations can be treated in the framework of 
the Hartree-Fock-Bogoliubov (HFB) approach and linear response theory. 

The neutron-drip line put a limit to the neutron-rich nuclei which can be 
produced in the laboratory or in supernova explosions (via rapid neutron- 
cupture processes). However, this limit can be overpassed in the inner crust 
of neutron stars since there the driped neutrons are kept together with the 
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neutron-rich nuclei by the gravitational preasure. 

The superfluid properties of the inner crust have been considered long 
ago in connection with post-glitch timing observations and colling processes 

1 O Q 

' ' . However, although the neutron star matter superfluidity has been 
intensively studied in the last decades ^, so far only a few microscopic 
calculations have been done for the superfluidity of inner crust matter. The 
most sophisticated microscopic calculations done till now use the framework 
of HFB approach at finite temperature ^. They will be discussed in the 
second part of my lecture. The discussion will be focused on the effects 
induced by the pairing correlations on the specific heat and on the cooling 
time of the inner curst of neutron stars. 



2. Pairing correlations in exotic nuclei 
2.1. Continuum- HFB approach 

The tool commonly used for treating pairing correlations in exotic nuclei 
close to the drip lines is the Hartree-Fock-Bogoliubov (HFB) approach ^ . 
The novel feature of interest here is how within this approach one can treat 
properly the quasiparticle states belonging to the continuum spectrum . 
This issue is discussed below. 

The HFB equations for local fields and spherical symmetry have the 
form: 

h(r)-\ A(r) \ fUi(r)\ =E (U i (r)\ (J) 



A(r) -h(r) + \J \Vi(r) J V^W 

where A is the chemical potential, h(r) is the mean field hamiltonian and 
A(r) is the pairing field. The fields depend on particle density p(r) and 
pairing density n(r) given by: 

P(r) = ^ £(2^ + 1)^(^(7-), (2) 

i 

KW = ^£(2* + w(r)vs(r). ( 3 ) 

i 

In the calculations presented here the mean field is described with a Skyrme 
force and for the pairing interaction it is used a density- dependent force 
of zero range of the following form ^: 

V(v - r') =V [1- v(-) a ]S(r - r') = V eff (p(r))S(v - r>). (4) 
Pa 
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For this force the pairing field is given by A(r) = ^^-n(r). 

The HFB equations have two kind of solutions. Thus, between and —A 
the quasiparticle spectrum is discrete and both upper and lower components 
of the radial HFB wave function decay exponentially at infinity. On the 
other hand, for E > —A the spectrum is continuous and the solutions are: 

uij(E,r) = C[cos(8i J )ji{air) - sm(5y)nj(air)] , 

v lj (E,r) = D 1 h ( i +) (i(3 1 r) , (5) 

where j; and ni are spherical Bessel and Neumann functions repectively 
and Sij is the phase shift corresponding to the angular momentum (Ij). 
The phase shift is found by matching the asymptotic form of the wave 
function written above with the inner radial wave function. From the en- 
ergy dependence of the phase shift one can determine the energy regions 
of quasiparticle resonant states. In HFB they are of two types. A first type 
corresponds to the single-particle resonances of the mean field. A second 
kind of resonant states is specific to the HFB approach and corresponds 
to the bound single-particle states which in the absence of pairing corre- 
lations have an energy e < 2 A. Among all possible resonances, of physical 
interest are the ones close to the continuum threshold. An cxcmplc of such 
resonances is shown in section 2.4. 



2.2. Resonant states in the BCS approach 

The low-lying quasiparticle resonances can be also calculated in the frame- 
work of the BCS approach. The BCS approximation is obtained by neglect- 
ing in the HFB equations the non-diagonal matrix elements of the pairing 
field. This means that in the BCS limit one neglects the pairing correlations 
induced by the pairs formed in states which are not time-reversed partners. 

The extension of BCS equations for taking into account the effect of 
resonant states was proposed in Refs.^' ^ For the case of a general pairing 
interaction the corresponding resonant-BCS equations read ^: 

A * = E V iifj u i v i + E V ii,u^ueu I 9u{e)u u {e)v v {e)de , (6) 

A " = E Ke^JpM + E V ^ye„— I 9,' (e'K' W)v* (e')de' , 

(7) 
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N = J2 V "+J2 I 9 v {e)vl{e)de . (8) 

i v 

Here Aj is the gap for the bound state i and A u is the averaged gap for 
the resonant state v. The quantity g v {e) = 2j, ^ +1 ^ is the continuum level 
density and 5 V is the phase shift of angular momentum {l v j u ). The factor 
g v {e) takes into account the variation of the localisation of scattering states 
in the energy region of a resonance ( i.e., the width effect) and becomes 
a delta function in the limit of a very narrow width. In these equations 
the interaction matrix elements are calculated with the scattering wave 
functions at resonance energies and normalised inside the volume where 
the pairing interaction is active. 

The BCS equations written above have been solved with a single particle 
spectrum corresponding to a HF ^ and a RMF ^ mean field. It was shown 
that by including in the BCS equations a few relevant resonances close 
to the continuum threshold one can get results very similar to the ones 
obtained in the HFB and RHB calculations. One can thus conclude that in 
nuclei close to the dripline the quasiparticle spectrum is dominated by a few 
low-lying resonances. In many odd-even nuclei close to the drip lines these 
low-lying quasiparticle resonances might be the only measurable excited 
states. Their widths can be calculated from the phase shift behaviour or 

19 

from the imaginary part of the energies associated to the Gamow states . 

In even-even nuclei the low-lying quasiparticle resonances could form 
the main component of unbound collective excitations with a finite time 
life. How these excitations can be treated in nuclei close to the drip line is 
discussed in the next section. 

2.3. Linear response theory with pair correlations and 
continuum coupling 

The collective excitations of atomic nuclei in the presence of pairing cor- 
relations is usually described in the Quasiparticle-Random Phase Approx- 
imation (QRPA) 6. In nuclei characterized by a small nucleon separation 
energy, the excited states are strongly influenced by the coupling with the 
quasiparticle continuum configurations. Among the configurations of par- 
ticular interest are the two-quasiparticle states in which one or both quasi- 
particlcs arc in the continuum. In order to describe such excited states one 
needs a proper treatment of the continuum coupling, which is missing in 
the usual QRPA calculations based on a discrete quasiparticle spectrum. 
In this section we discuss how the pairing and the continuum coupling can 
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be treated in the framework of the linear response theory 6 ' . 

The response of the nuclear system to an external perturbation is ob- 
tained from the time-dependent HFB equations (TDHFB) ^: 

ih ^ = m) + m ,n{t)], (9) 

where 1Z and Tt are the time-dependent generalized density and the HFB 
Hamiltonian, respectively. T is the external oscillating field : 

T = Fer iut + h.c. (10) 

In Eq. (10) F includes both particle-hole and two-particle transfer operators 



f = Y. 4Mc, + £(4 2c H + if CiCj ), (ii) 

ij ij 

and c\, Ci are the particle creation and annihilation operators, respectively. 
In the small amplitude limit the TDHFB equations become: 

hull' = [H\ Tl ] + [H°, W] + [F, 11% (12) 

where the superscript ' stands for the corresponding perturbed quantity. 

The variation of the generalized density W is expressed in term of 3 
quantities, namely p' , n' and R' , which are written as a column vector : 




where p'^ = ^0|cjci|'y is the variation of the particle density, n'^ = (0\cjCi\') 

and R'^ = (o\c^c\\'^ are the fluctuations of the pairing tensor associated 

to the pairing vibrations and |') denotes the change of the ground state 
wavefunction |0 > due to the external field. 

The variation of the HFB Hamiltonian is given by: 



n' = vp', (w) 

where V is the matrix of the residual interaction expressed in terms of the 
second derivatives of the HFB energy functional, namely: 

V^(ra,rV)= a, (3 =1,2,3. (15) 
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In the above equation the notation a means that whenever a is 2 or 3 then 
a is 3 or 2. 

Introducing for the external field the three dimensional column vector: 




(16) 

the density changes can be written in the standard form: 

p' — GF , (17) 

where G is the QRPA Green's function obeying the Bethe-Salpeter equa- 
tion: 



G = (1 - Gov)" 1 G = G + G VG. (18) 
The unperturbed Green's function Go has the form: 



G QP (r<7,iV;w) 



U^{va)U^{v'a') ^ 2 (ra)^ 2 (rV')_ 

■V 

(19) 



Kuj - (Ei +Ej) + iT] hjj + (Ei +E j ) + iTj 

(1 

where Ei are the qp energies and Uij are 3 by 2 matrices expressed in term 
of the two components of the HFB wave functions The J] symbol in 
Eq. (19) indicates that the summation is taken over the discrete and the 
continuum quasiparticle states. 

The QRPA Green's function can be used for calculating the strength 
function associated with various external perturbations. For instance, the 
transitions from the ground state to the excited states induced by a particle- 
hole external field can be described by the strength function: 



S(w) = -^Im J F 11 *(r)G 11 (r,r / ;w)J ,11 (r / )dr dr' (20) 

where G 11 is the (ph,ph) component of the QRPA Green's function. An- 
other process which can be described in the same manner is the two-particle 
transfer from the ground state of a nucleus with A nucleons to the excited 
states of a nucleus with A+2 nucleons. For this process the strength func- 
tion is: 
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S(u) = ~Im J F 12 *(r)G 22 (r,r';w)F 12 (r')dr fir' (21) 
where G 22 denotes the (pp,pp) component of the Green's function. 

2.4. Quasiparticle excitations and two-neutron transfer in 
neutron-rich nuclei 

The formalisms presented above are illustrated here for the case of neutron- 
rich oxygen isotopes First, we present an exemple of quasiparticle 
resonances calculated in the framework of the continuum-HFB (cHFB) ap- 
proach introduced in section 3.1. In the cHFB calculations the mean field 
quantities are evaluated using the Skyrme interaction SLy4 ^ , while for 
the pairing interaction we take a zero-range density-dependent force. The 
parameters of the pairing force are given in Ref The HF single-particle 
and HFB quasiparticle energies corresponding to the sd shell and to the 
I/7/2 state are listed in Table 1. One can notice that in both HF and cHFB 
calculations the state I/7/2 is a wide resonance for 18_22 nuclei, while the 
state ld 3 / 2 is a narrow resonance. As seen below, these one-quasiparticle 
resonances determine essentially the two-quasiparticle states which are the 
stongest populated in even-even oxygen isotopes. 



Table 1. HF and HFB energies in oxygen isotopes. 







18Q 


20Q 


22Q 


ld 5 /2 


HF 


-6.7 


-7.0 


-7.45 


ld 5 /2 


cHFB 


2.26 


2.08 


2.30 


2si/ 2 


HF 


-4.0 


-4.2 


-4.6 


2Sl/2 


cHFB 


3.46 


2.28 


1.05 


ld 3 /2 


HF 


(0.46;0.02) 


(0.51;0.03) 


(0.42;0.02) 


ld 3/2 


cHFB 


(7.74;0.12) 


(6.60;0.29) 


(5.39;0.01) 


lf 7/ 2 


HF 


(5.50;1.35) 


(5.24;1.24) 


(4.86;1.04) 


1I7/2 


cHFB 


(12.77;1.13) 


(12.14;0.83) 


(10.05;0.69) 



The two-quasiparticle states are calculated by using the response theory 
described in section 2.3. In the calculations one includes the full discrete 
and continuum HFB spectrum up to 50 MeV. These states are used to 
construct the unperturbed Green's function Go- After solving the Bcthc- 
Salpeter equation for the QRPA Green function one constructs the strength 
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functions written in section 2.3. The details of the calculations can be found 
m Ref 14 . Here we show the results obtained for the collective states excited 
in the two-neutron transfer. The strength function corresponding to a neu- 
tron pair transferred to the oxygen isotope 22 O is shown in Fig.l. 




Fig. 1. The response function for the two-neutron transfer on 220. The unperturbed 
response is in solid line and the QRPA response in dashed line. 



For the isotope 22 O the subshell c? 5 / 2 is essentially blocked for the pair 
transfer. Therefore in this nucleus we can clearly identify only two peaks 
below 11 MeV, corresponding to a pair transferred to the states 2s\/2 and 
2c?3/2- The strength function shown in Fig. 1 shows also a broad peak around 
20 MeV. This peak is built mainly upon the single-particle resonance I/7/2 
and its cross section is much larger than the one associated to the lower en- 
ergy transfer modes. Since this high energy transfer mode is formed mainly 
by single-particle states above the valence shell, this mode is similar to 
the giant pairing vibration mode suggested long ago Although such a 

mode has not been detected yet, the pair transfer reactions involving exotic 

1 (\ 

loosely bound nuclei may offer a better chance for this undertaking . 
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3. Superfluid and thermal properties of neutron stars crust 

In this section we analyse the superfluid properties of a nuclear system in 
which the limit of the neutron-drip line is overpassed: the inner crust of 
neutron stars. The inner crust consists of a lattice of neutron-rich nuclei 
immersed in a sea of unbound neutrons and relativistic electrons. Down to 
the inner edge of the crust, the crystal lattice is most probably formed by 
spherical nuclei. More inside the star, before the nuclei dissolve completely 
into the liquid of the core, the nuclear matter can develop other exotic 
configurations as well, i.c, rods, plates, tubes, and bubbles °. The thickness 
of the inner crust is rather small, of the order of one kilometer, and its mass 
is only about 1% of the neutron star mass. However, in spite of its small 
size, the properties of inner crust matter, especially its superfluid properties, 
have important consequences for the dynamics and the thermodynamics of 
neutron stars. In what follows we will discuss the main features of the 
pairing correlations inside the inner crust matter. Then we shall focus on 
the effects induced by the nuclear superfluidity on the specific heat and on 
the cooling time of the curst. 



3.1. Superfluid properties of the inner crust matter 

The superfluid properties of the inner crust matter discussed here are 
based on the finite-temperature HFB (FT-HFB) approach. For zero range 
forces and spherical symmetry the radial FT-HFB have a similar form as 
the HFB equations at zero temperature, ie., 

h T (r)-X A T (r) \[U i (r)\_ [U i (r)\ 
A T (r) -h T (r)+\)\Vdr)) ^ \V % {r) ) ' ^ 

where E{ is the quasiparticle energy, Ui , Vi are the components of the radial 
FT-HFB wave function and A is the chemical potential. The quantity h,T(r) 
is the thermal averaged mean field hamiltonian and Ax(r) is the thermal 
averaged pairing field. The latter depends on the average pairing density 
kt- In a self-consistent calculation based on a Skyrme-type force , hrir) is 
expressed in terms of thermal averaged densities, i.e., kinetic energy density 
tt(t), particle density prir) and spin density Jr(r), in the same way as in 
the Skyrme-HF approach. The thermal averaged densities mentioned above 
are given by . 

Mr) = ^ + l)M*(rW(r)(l - .A) + U* (r^r) fr] 
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^(r) = ± £(2* + l)i;;(r)^(r)(l - 2/,) 

i 

Mr) = ^ + Vl^Ji + 1) - Hh + 1) - J] 

i 

x{^ 2 (l-/i) + f/f/i} 

i 

x(l-/.) + l(^-?)' + ^0fl/.). 

(23) 

where /« = [1 + ea^-Ei/fcBT)] -1 is the Fermi distribution, fcs is the Boltz- 
mann constant and T is the temperature. The summations in the equations 
above are over the whole quasiparticle spectrum, including the unbound 
states. 

The FT-HFB equations are solved for the spherical Wigner-Seitz cells 
on 

determined in Ref. . To generate far from the nucleus a constant density 
corresponding to the neutron gas, the FT-HFB equations are solved by 
imposing Dirichlct-Ncumann boundary conditions at the edge of the cell 

on 

, i.e., all wave functions of even parity vanish and the derivatives of odd- 
parity wave functions vanish. In the FT-HFB calculations we use for the 

1 7 

particle-hole channel the Skyrmc effective interaction SLy4 ' , which has 
been adjusted to describe properly the mean field properties of neutron- 
rich nuclei and infinite neutron matter. In the particle-particle channel we 
employ a density dependent zero range force. Since the magnitude of pairing 
correlations in neutron matter is still a subject of debate, the parameter 
of the pairing force are chosen so as to describe two different scenarios for 
the neutron matter superfluidity Thus, for the first calculation we use the 
parameters: Vb=-430 McV fm 3 , 77=0.7, and a=0.45. With these parameters 
and with a cut-off energy for the quasiparticle spectrum equal to 60 MeV 
one obtains approximately the pairing gap given by the Gogny force in 
nuclear matter ^. In the second calculation we reduce the strength of the 
force to the value Vb=-330 MeV fm 3 . With this value of the strength we 
simulate the second scenario for the neutron matter superfluidity, in which 
the screening effects would reduce the maximum gap in neutron matter to 
a value around 1 MeV . 

The FT-HFB results are shown here for two representative Wigncr- 
Seitz cells chosen from Ref. . These cells contain Z=50 protons and have 



February 9, 2008 0:50 



WSPC - Proceedings Trim Size: 9in x 6in nsandulescu 



11 



rather different baryonic densities, i.e., 0.0204 fm -3 and 0.00373 fm -3 . The 
cells, which contain N=1750 and N=900 neutrons, respectively, are denoted 
below as a nucleus with Z protons and N neutrons, i.e., 1800 Sn and 950 Sn. 
The FT-HFB calculations are done up to a maximum temperature of T=0.5 

99 

MeV, which is covering the temperature range of physical interest . 




Fig. 2. Neutron pairing fields for the cell 800 Sn calculated at various temperatures. 
The numbers 1 and 2 which follow the cell symbol (see the inset) indicate the vari- 
ant of the pairing force used in the calculations. The full and the dashed lines cor- 
responds (from bottom upwards) to the set of temperatures T={0.0, 0.5} MeV and 
T={0.0, 0.1, 0.3, 0.5}MeV, respectively. 



The temperature dependence of the pairing fields in the two cells pre- 
sented above is shown in Figs. 2-3. First, we notice that for all temperatures 
the nuclear clusters modify significantly the profile of the pairing field. One 
can also see that for most of the cases the temperature dependence of the 
pairing field is significant. This is clearly seen in the low-density cell 950 Sn. 

3.2. Specific heat of the inner crust baryonic matter 

The supcrfluid properties of the neutrons discussed in the previous section 
have a strong influence on the specific heat of the inner crust matter ^ . The 
specific heat of a given cell of volume V is defined by: 



1 d£{T) 
~ V dT 



(24) 
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Fig. 3. The same as in Fig.l, but for the cell 950 Sn. The full and the dashed lines 
corresponds (from bottom upwards) to the set of temperatures T={0.0, 0.1, 0.3, 0.5} and 
T={0.0,0.1}MeV, respectively 



where £ (T) is the total energy of the baryonic matter inside the cell, i.e., 



Due to the energy gap in the excitation spectrum, the specific heat of 
a superfluid system is dramatically reduced compared to its value in the 
normal phase. Since the specific heat depends exponentially on the energy 
gap, its value for a Wigner-Seitz cell is very sensitive to the local variations 
of the pairing field induced by the nuclear clusters. This can be clearly 
seen in Fig. 4, where the specific heat is plotted for the cell 1800 Sn and for 
the neutrons uniformly distributed in the same cell. One can notice that 
at T=0.1 MeV and for the first pairing force the presence of the cluster 
increases the specific heat by about 6 times compared to the value for the 
uniform neutron gas. However, the most striking fact seen in Fig. 4 is the 
huge difference between the predictions of the two pairing forces. Thus, for 
T=0.1 MeV this difference amounts to about 7 orders of magnitude. 

The behaviour of the specific heat for the low- density cell 950 Sn is shown 
in Fig. 5. For the first pairing force we can also see that at T=0.1 MeV the 
cluster increases the specific heat by about the same factor as in the cell 
1800 Sn. However, for the second pairing force the situation is opposite: the 
presence of the nucleus decreases the specific heat instead of increasing it. 




(25) 
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Fig. 4. Specific heat for the cell 1800 Sn as a function of temperature. The notations used 
in the inset and the representation of the calculated values are the same as in Figs. 1-3. 




-16 

0.1 0.2 0.3 

T [MeV] 

Fig. 5. The same as in Fig. 4, but for the cell 950 Sn. 



3.3. Collective modes in the inner crust matter 

In the calculations presented in the previous section the specific heat of in- 
ner crust matter was evaluated by considering only non-interacting quasi- 
particles states. However, the specific heat can be also strongly affected 
by the collective modes created by the residual interaction between the 
quasiparticles, especially if these modes appear at low-excitation energies. 
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The collective modes in the inner crust matter were calculated in Ref. 

in the framework of linear response theory discussed in section 2.3. The 
most important result of these calculations is the apparence of very col- 
lective modes at low energies, of the order of the pairing gap. An example 
of such mode is seen in Fig. 6, were is shown the quadrupole response for 
the cell 1800 Sn. As can be clearly seen from Fig. 6 , when the residual inter- 
action is introduced among the quasiparticles the unperturbed spectrum, 
distributed over a large energy region, is gathered almost entirely in the 
peak located at about 3 MeV. This peak collects more than 99% of the 
total quadrupole strength and is extremely collective. An indication of the 
extreme collectivity of this low-energy mode can be also seen from its re- 
duced transition probability, B(E2), which is equal to 25. 10 3 Weisskopf 
units. This value of B(E2) is two orders of magnitude higher than in stan- 
dard nuclei. This underlines the fact that this WS cell cannot be simply 
considered as a giant nucleus. The reason is that in this cell the collective 
dynamics of the neutron gas dominates over the cluster contribution. 
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Fig. 6. Quadrupole strength distribution of neutrons for the cell Sn. The full curve 
represents the QRPA strength, and the dashed line is the HFB unperturbed strength. 
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The collective excitations located at low energies can affect significantly 
the specific heat of the inner crust baryonic matter. This can be seen in 
Fig. 7, where the specific heat corresponding to the collective modes (of mul- 
tipolarity L=0,l,2,3) is shown. We notice that for T=0.1 MeV the specific 
heat given by the collective modes is of the same order of magnitute as 
the one corresponding to HFB spectrum. Therefore one expects that the 
collective modes of the inner crust matter could affect strongly the thermal 
behaviour of the crust. 



-13 




0.1 0.2 0.3 0.4 0.5 

T [MeV] 

Fig. 7. Specific heat in the cell 1800 Sn. The dashed line corresponds to the collective 
modes and the full line to the HFB spectrum 



3.4. Cooling time of the inner crust 

The specific heat of the inner crust matter is an important quantity for 
cooling time calculations. In this section we shall discuss the sensivity of 
the cooling time on the specific heat calculated with the two scenarios for 
the nuclear superfluidity introduced in section 3.1. The cooling process we 
analyse here corresponds to a fast cooling mechanism (e.g., induced by 
direct Urea reactions) . In this case the interior of the star cools much faster 
than the crust. The cooling time is defined as the time needed for the cooling 
wave to traverse the crust and to arrive at the surface of the star. According 
to numerical simulations , the cooling time is proportional to the square 
of the crust size and to the specific heat. This result was used by Pizzochero 
ct al ^ for estimating the cooling time in a simple model, which we have 
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also employed in our calculations . Thus, the crust is devided in shells 
of thickness Ri for which the thermal difusivity Di could be considered as 
constant. Then the cooling time is obtained by summing the contribution 
of each shell, i.e., 



The thermal diffusivity is given by D = where n is the thermal con- 
ductibility and Cy is the specific heat. The thermal conductivity is mainly 
determined by the electrons and in our calculations we have used the values 
reported by Lattimer et al. . The specific heat Cy has major contributions 
from the electrons, which can be easily calculated, and from the neutrons 
of the inner crust. As we have seen above, the specific heat of the neutrons 
depends strongly on pairing correlations. In order to see if this dependence 
has observationally consequences for the cooling time, we have performed 
two calculations, corresponding to the strong and the weak pairing forces 
introduced in section 3.1. The calculations show that the cooling time is 
increasing with more than 80% if for calculating the pairing correlations 
we use a weak pairing force instead of a strong pairing force. This result 
indicates that the nuclear superfluidity of the inner crust matter plays a 
crucial role for the cooling time calculation of neutron stars. 

4. Conclusions 

In this lecture we have shown how the HFB approach and the linear re- 
sponse theory can be used to describe the pairing correlations in exotic 
nuclei and in the inner crust of neutron stars. Thus, for the nuclei close 
to the drip lines we have discussed how one can incorporate in the two 
approaches mentioned above the effects of the continuum coupling on pair- 
ing correlations. Then, using the same models, we have analysed what are 
the effects of pairing correlations on the specific heat and on the cooling 
time of neutron stars. We thus found that the cooling time depends very 
strongly on the intensity of pairing correlations in nuclear matter. Because 
the intensity of pairing correlations in nuclear matter is still unclear, at 
present one can only estimate the limits in which the cooling time of the 
inner crust can vary. Since the pairing correlations in nuclear matter and in 
nuclei are in fact correlated, one hopes that these limits could be reduced 
by systematic studies of pairing in both infinite and finite nuclear systems. 
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